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BABAR was the first large experiment to incorporate Geant4 into its detector simulation. Since July 2001, 
1.5x 10 9 BABAR events have been produced using this simulation. In a typical e + e — — > T(4s) — > B°B° event, 
between 30 and 60 tracks are produced in the generator and propagated through the detector, using decay, 
electromagnetic and hadronic processes provided by the Geant4 toolkit. The material model of the detector is 
very detailed and a special particle transportation module was developed so that minute features (on the few 
micron scale) would be sampled in the propagation without sacrificing performance. The propagation phase 
for such an event requires 5 CPU seconds on an 866 MHz processor. Execution speeds for other BABAR event 
types will also be presented. Validation of simulated events with BABAR data is ongoing, and results of Monte 
Carlo/data comparisons will be shown. A discussion of the design of the simulation code, how the Geant4 toolkit 
is used, and ongoing efforts to improve the agreement between data and Monte Carlo will also be presented. 



1. OVERVIEW OF THE BaBar 
SIMULATION 

1.1. The Detector 

The primary physics goal of the BABAR experi- 
ments the study of CP violation in the B°B° system. 
B°W events are produced by the decay of T(4s) res- 
onances which are produced in e + e~ collisions. Some 
of the B°s decay into CP eigenstates such as J/'fAT" , 
and the subsequent decays of the J/^f and are 
detected and reconstructed. A typical event will pro- 
duce between 30 and 60 tracks, and typical energies 
of the final state decay products are: 

• lepton pairs: 0.3 < p < 2.3 GeV/c 

• tt° : 0.3 < E < 2.5 GeV 

• 7 : 0.1 < E < 4.5 GeV . 

Hadronic final states are also important because 
charged 7rs and Ks interact in the beam pipe and 
calorimeters. Typical hadron momenta are 

• p < 4 GeV/c, with most < 1 GeV/c . 

Background events such as Bhabha scattering produce 
the highest energy tracks, with p e - ,p e + < 9 GeV/c . 

The BABAR detector was built to collect final state 
tracks with high efficiency and precision, and to al- 
low the reconstruction of B° decays into a wide range 
of exclusive final states with low background. The 
detector was also designed to operate at the high lu- 
minosities provided by the PEP-II accelerator. 



The detector consists of a silicon vertex tracker 
(SVT) which surrounds the e + e~ interaction point 
and provides vertex determination at the 10 /im level. 
Surrounding the SVT is a He-isobutane drift cham- 
ber (DCH) for the measurement of charged tracks. 
Outside the DCH is a detector of internally-reflected 
Cherenkov radiation (DRC), constructed of quartz 
bars and used for particle identification (PID). Out- 
side the DRC is a Csl electromagnetic calorimeter 
(EMC) designed to contain showers from charged and 
neutral tracks. A muon tracker and instrumented flux 
return (IFR), built of alternating layers of iron and re- 
sistive plate chamber (RPC) detectors, surrounds all 
the interior detectors. 

1.2. The Simulation 

The simulation was designed to satisfy four main 
requirements. It must: 

• run within the BABAR Framework 0. All the 
tracking, physics and hit-scoring performed by 
Geant4 |2| is implemented as a Framework mod- 
ule. Because the Framework is responsible for 
run and event control, Geant4 must therefore 
relinquish this duty. 

• work with existing event generators, detector re- 
sponse codes and reconstruction codes. At the 
time Geant4 became available, BABAR detector 
response and reconstruction codes were essen- 
tially complete and many event generation codes 
were already available. 
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Figure 1: Simple diagram of the BABAR simulation. Boxes above the horizontal line represent BABAR Framework 
modules, and the wide arrows depict the event flow through the Framework. Narrow arrows indicate the retrieval of 
information from the database. 



use the Objectivity database Q for persistence. 2. VALIDATION 



• be detailed, yet fast enough to keep up with 
high-luminosity event production. 

A simple diagram of the main functions and opera- 
tion of the resulting simulation is shown in Fig. ^ 

The Geant4 toolkit is used extensively in the simu- 
lation. The detector geometry is built with Geant4 
simple volumes and Boolean operations. As men- 
tioned above, Geant4 is also responsible for hit-scoring 
and provides all the physics processes. These include 
the decay, standard electromagnetic and low-energy 
hadronic [E < lOGeV) processes. Some features pro- 
vided by Geant4, such as detector response code and 
persistence, are not used. The default Geant4 trans- 
portation/navigation process is also not used. It em- 
ploys a general Runge-Kutta stepper which was found 
to be too slow and not precise enough for the many 
thin volumes in the BABAR simulation. Taking ad- 
vantage of the slowly varying magnetic field in the 
BABAR detector, a specialized stepper was developed 
to meet these needs. This stepper determines how far 
a particle can travel before the field deviates from a 
locally constant value. Perfect helices are then used 
to propagate the particle over that distance. If that 
distance is large, many small steps are avoided. 



Before being put into production, the 
BABARGeant4 simulation was subjected to a 
series of validation tests focusing on 

• verifying the detector material model, 

• electromagnetic processes, 

• hadronic processes, 

• tracking, resolution and reconstruction, 

• particle ID, 

• performance and robustness. 

Since October 2000, 20 million simulated events of 
several types have been produced for these tests. They 
include generic B°B°, Bhabha scattering and dimuon 
events, among others. Comparsion of these events 
with data has allowed a refinement of the detector 
material model to the point that no further changes 
are necessary. Understanding of particle identifica- 
tion and electromagnetic processes is also believed to 
be well in hand, although validation in these areas 
continues. However, additional validation is required 
for the hadronic processes. 



TUMT006 



Computing in High Energy Physics, La Jolla, CA, March 24-28, 2003 



SLAC-PUB-9862 3 




ADC counts 



850 900 



9 



Figure 2: dE/dx for minimum ionizing electrons and positrons from Bhabha scattering. The histograms represent two 
different simulations of the energy loss in the drift chamber gas. The points are from a BABAR data run. The 
histograms are normalized to the number of events in the data run. 



2.1. Electromagnetic Process Validation 

The BABAR simulation uses the following Geant4 
electromagnetic processes 

• for photons: photo-electric effect, Compton 
scattering and e + e~ pair production; 

• for electrons and positrons: electron ionization, 
electron bremsstrahlung and e + annihilation; 

• for muons: [i ionization, bremsstrahlung and 
e + e~ pair production by muons; 

• for hadrons: hadron ionization; 

• for all charged particles: multiple scattering. 

One of the most basic validation tests is to repro- 
duce the expected energy loss due to ionization pro- 
cesses. In this case, minimum ionizing electrons and 
positrons from radiative Bhabha scattering were ex- 
amined in the He-Isobutane gas of the drift cham- 
ber. Simulated events are compared to data from a 
recent run in Fig. [21 where it is seen that both the 
shape and normalization of the dE/dx distributions 
are in agreement to within 15%. This indicates that 



fluctuations in the energy loss are reasonably well- 
reproduced. More importantly, the mean values of 
dE/dx agree to within 1%. 

It is also important to study electromagnetic show- 
ers in the Csl calorimeter. Measuring shower shapes 
provides a cumulative test of the photon and electron 
processes listed above. Radiative Bhabha scattering 
was used to examine the distribution of azimuthal 
and lateral shower parameters as a function of elec- 
tron/positron energy. Fig. |31 shows the lateral distri- 
bution parameter (LAT). With the exception of a few 
energy bins, distributions of both of these parameters 
were found to show good agreement between data and 
simulation for all electron and positron energies up to 
4 GeV. 

2.2. Reconstruction 

The comparison of reconstructed events from data 
and simulation provides a stringent validation test. 
In BABAR an accurate reconstruction of the 7r° mass 
and width provides a test of many detector per- 
formance features, including tracking, energy scale, 
shower development, shower containment and detec- 
tor response. 
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Figure 3: Distribution of the LAT parameter as a function of incident electron or positron energy. Top left: all 
energies; top right: to 1 GeV; bottom left: 1 to 3 GeV; bottom right: 3 to 5 GeV. The histograms represent a 
simulated sample of Bhabha events and the points are taken from BABAR data. 
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Figure 4: 7r° mass reconstruction for data (left) and simulation (right). Histograms represent the reconstructed samples 
and the curves are fits to the samples. The total pion energy is between 0.8 and 1.0 GeV. The fitted mass and width 
values are stable over the range 0.3 - 2.1 GeV. 



In B°B° events, the decay K s — > 7r°7r° provides 
7r°s with energies between 0.3 and 2.1 GeV. Fig. 0] 
shows the reconstructed mass for ir°s with total ener- 
gies between 0.8 and 1.0 GeV. For both data and sim- 
ulated samples the mass is close to the expected value, 



which indicates that tracking, energy scale and shower 
containment are well-understood. The widths, how- 
ever, are significantly different and point to problems 
in shower development or detector response. Since 
it is believed that the shower development is under- 
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Figure 5: Cross section for ir + production from 730 MeV protons on carbon at 15°, 45°, 90° and 135°. Histogram: 
prediction of the Low Energy Parameterized (LEP) model in Geant4. Data points are from Ref. 



stood, detector response should be examined in future 
validation tests. 



2.3. Hadronic Process Validation 

The electromagnetic processes of hadron ionization 
and multiple scattering are of primary importance in 
the propagation of long-lived hadrons through the de- 
tector. Hadronic interactions are still important, how- 
ever, because scattering from nuclei and the produc- 
tion of hadronic secondaries in the dense material of 
the EMC and IFR can affect the energy deposition. 

Hadronic interactions in the BABAR simulation are 
currently handled by the so-called "low energy pa- 
rameterized" , or LEP, model. This is a re-engineered 
version of the GHEISHA code 0, which has been 
tuned for hadrons with incident energies below 20 



GeV. Even so, it is not particularly appropriate at 
BABAR hadron energies which are typically below 1 
GeV. This is demonstrated in Fig. |5] which shows the 
pion production cross section for 730 MeV protons on 
a thin carbon target. Large, qualitative differences be- 
tween the LEP model and data are evident. Recently, 
however, better models, such as the Bertini cascade 
0, have become available. As shown in Fig. El much 
better agreement with the data is achieved at all an- 
gles. A current limitation of this model is that it is 
not valid for kaons. 

The hadronic process validation in BABAR has so 
far concentrated on the comparison of results from 
various models to published data from other exper- 
iments in which the incident hadron energies are in 
the BABAR range. Validations using BABAR data are 
just beginning. 
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Figure 6: Cross section for ir + production from 730 MeV protons on carbon at 15°, 45°, 90° and 135°. Histogram: 
prediction of the Bertini cascade model in Geant4. Data points are from Ref. |5J. 



While BABAR is certainly not optimized for 
hadronic physics validation, two regions of the detec- 
tor are useful for tests. The beam pipe support tube 
is a cylinder of carbon fiber and the inner wall of the 
drift chamber is a cylinder of beryllium. Tests are just 
beginning which use these regions as thin targets for 
incident hadrons. 



Table I Execution Time for BABAR Events 



Event type 


CPU time (sec) 


B0B0 


5.0 


bhabha 


7.0 


dimuon 


0.6 



3. PERFORMANCE 

The simulation stage of BABAR event production 
includes event generation, tracking and hit-scoring. A 
Pentium III 866 MHz PC currently requires 5.0 sec- 
onds on average to simulate a B°B° event. Average 
execution times for other BABAR event types are given 
in Table IJ 



Up to now, relatively little effort has been devoted 
to optimizing performance because validation was the 
main concern. However, some increase in speed is 
expected with 

• improvements in BABAR geometry models, 

• improvements in the Geant4 electromagnetic 
processes, and 
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• fine tuning of the secondary particle production 
cuts. 



overall speed of the simulation. 



On the other hand, the improved hadronic models 
are much more CPU-intensive than the currently used 
LEP model, and will cause some slowdown. 

After the production of 20 million validation events, 
many of the bugs in both the BABAR simulation code 
and in Geant4 have been shaken out. To date, BABAR 
simulation production has generated more than 1.5 
billion events of all types at roughly 20 production 
sites in North America and Europe, with a failure rate 
of less than one event per million in the simulation 
stage. 



4. CONCLUSION 

BABAR is the first large experiment to develop and 
use a Geant4-based simulation. A detailed detector 
model has been developed using the Geant4 geome- 
try, and both electromagnetic and hadronic interac- 
tions are implemented with Geant4 physics processes. 
The simulation was designed to run within the BABAR 
Framework application. 

Comparisons of simulated event samples and 
BABAR data indicate that the electromagnetic pro- 
cesses perform largely as expected in the few GeV 
energy range. Validation tests of the LEP hadronic 
model show that more detailed models will be required 
before good agreement with data can be achieved. 

The BABAR simulation has proven to be robust, 
with 1.5 billion events generated so far at a low failure 
rate. Future efforts will concentrate on improving the 
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